Skyrmion Physics Beyond the Lowest Landau Level Approximation 
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The effects of Landau level mixing and finite thickness of the two-dimensional electron gas on the 
relative stability of skyrmion and single spin-flip excitations at Landau level filling factor v = 1 have 
been investigated. Landau level mixing is studied by fixed-phase diffusion Monte Carlo and finite 
thickness is included by modifying the effective Coulomb interaction. Both Landau level mixing 
and finite thickness lower skyrmion excitation energies and favor skyrmions with fewer spin flips. 
However, the two effects do not work 'coherently'. When finite thickness is included the effect of 
Landau level mixing is strongly suppressed. 
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In a seminal paper Sondhi et al.a predicted that if the 
ratio of the Zeeman energy to the exchange energy is suf- 
ficiently small the quasiparticle excitations of the ly = 1 
integer quantum Hall state should be char ged-j skyrmions' 
with topologically nontrivial spin texturesu The exis- 
tence of these [exotic excitations is now well confirmed 
by experiment .□□ 

While the initial theoretical description of skyrmions 
in terms of an effective long-wavelength nonlinear sigma 
model is appropriate for excitations whose magnetiza- 
tion varies slowly on the scale of the magnetic length, 
in typical experiments the number of reversed spins per 
skyrmion is ~ 3 — 4 and a more microscopic descrip- 
tion is required. A step towards such a description was 
taken by MacDonald, Fertig, and BreyH who introduced 
explicit wave functions describing skyrmions with well 
defined spin quantum numbers. 

In this paper we present the results of calculations in 
which these microscopic wave functions have been used to 
study the effects of finite thickness of the two-dimensional 
electron gas (2DEG) and Landau level mixing on ly — 1 
skyrmion states. Finite thickness is included by modi- 
fying the Coulomb interaction and Landau level mixing 
is studied-, using the fixed-phase diffusion Monte Carlo 
method.EfQ One motivatipn for this work is the observa- 
tion by Schmeller et al.u of skyrmion energy gaps mea- 
sured in transport experiments which are 25% smaller 
than theoretical estimates based on lowest Landau level 
calculations ,eI a result which led these authors to spec- 
ulate that Landau level mixing plays an important role 
in reducing skyrmion excitation energies. In our calcula- 
tions we find that both finite thickness and Landau level 
mixing do in fact lower skyrmion excitation energies and 
favor skyrmions with fewer spin flips. However, the two 
effects do not work efficiently together — the finite thick- 
ness effect strongly suppresses the Landau level mixing 
effect. 

In our, calculations we have used the spherical 
geometrya in which N electrons with charge — e and mass 
m are confined to move on the surface of a sphere of ra- 
dius R in the presence of a radial magnetic field B. The 
total number of flux quanta piercing the surface of the 



sphere is then N^ — 2 (R/Iq) , where ^o = \/hc/eB is 
the magnetic length, and the Hamiltonian describing the 
system is 
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Vc + Hz. (1) 



Here Vc = J2i<i ^i^ij) i^ the Coulomb interaction where 
rij is the chord distance between electrons i and j and 
Hz — —gfJ-B J2i SizB is the Zeeman interaction. We 
work in the Wu-Yang gauge in which the, vector poten- 
tial is A = e0?ic7V0(l-cos6l)/(2ei?sin6').N The softening 
of the short-range part of the Coulomb interaction due 
to finite thickness of the lowest subband wave function of 
the 2DEG is included through the modified interaction 
V{r) — e^/e-\/r^-|-pS? where /3 is a parameter character- 
izing the thickness .O 

For the v = 1 ground state the electron spins are 
fully polarized and the lowest Landau level is completely 
full. On the sphere this occurs when N = N^ + 1 and 
the corresponding wave function is a Slater determinant 
of lowest Landau level wave functions. If the electron 
positions are given in terms of the complex coordinate 
z = tan(0/2) exp (— 1(/)), where 9 and (f) are the usual 
spherical angles, then, exploiting the Vandermonde form 
of the Slater determinant, the ground state wave function 
can be written 
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where C/ = nfc(l + Nfe|')"'- 

Single spin-flip excitations of the system are con- 
structed by removing a spin up electron from the ground 
state, flipping its spin, and reintroducing it into the low- 
est spin-down Landau level. Letting K denote the num- 
ber of reversed spins associated with a given quasiparti- 
cle excitation this procedure produces a K = quasihole 
with charge +e and a K = 1 quasielectron with charge 
— e. The transport gap which determines the activated 
temperature dependence of the longitudinal resistivity at 
:^ = 1 is set by the excitation energy for creating a quasi- 
electron and quasihole with infinite separation. On the 



sphere the best approximation to this state can be real- 
ized by placing the quasihole at the bottom of the sphere 
{9 = tt) and the quasielectron at the top of the sphere 
{0 = Q). The corresponding wave function is 



V-,^ = A C/"*/2 
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where the operator A antisymmetrizes the wave function 
under the exchange of all pairs of electrons. 

The nature of the quasiparticle excitations at ly = 1 
depends on the relative size of the Zeeman energy and 
the exchange energy characterized by the dimensionless 
ratio g = gfj,BB/{e'^/elo). In the limit of large g the 
transport gap for the ly = 1 integer quantum Hall state 
is set by the single spin-flip excitation described above. 
However, as g is lowered below a critical value the single 
spin flip becomes unstable to the formation of a neutral 
skyrmion pair, i.e. a well separated charge +e and charge 
— e skyrmion pair, involving more than one flipped spin.EJ 
A microscopic wave function description of the charge +e 
skyrmion excitations was developed by MacDonald, Fer- 
tig, and BreyB who showed that for a model Hamiltonian 
with a hard-core delta-fmiction repulsion the states 

,PskiK)^AU^''^+''>^'l[{z,-z,) 
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are exact eigenstates of H with charge e and K reversed 
spins. 

The wave functions ijjgs, ipsf and ipsk{K) describe 
states in which the kinetic energy is completely quenched 
and there is no Landau level mixing. Therefore the only 
contributions to the energy gaps computed using these 
wave functions will come from the Coulomb and Zeeman 
energies. The variational result for the Coulomb contri- 
bution to the single spin-flip energy gap is 



A,/ = {^sf\Vc\A}) - {^gs\Vc\^g 



(5) 



In the absence of Landau level mixing an exact particle- 
hole symmetry in the lowest Landau level transforms a 
charge +e skyrmion with K reversed spins into a charge 
— e skyrmion with K -\- 1 reversed spins. It is therefore 
possible to compute the energy gap for creating a neu- 
tral skyrmion pair by first computing A^/ and then com- 
puting the change in Coulomb energy of the charge +e 
skyrmion as the number of reversed spins is increased, 

5{K) - {iysk{K)\Vc\i'sk{K)) - (V'sfe(O)ll^cl^sfc(O)). (6) 

As pointed out by Abolfath et al.E3 the wave functions 
(Q) involve power-law tails which lead to strong finite-size 
effects. In our variational calculations we have computed 
5{K) for /3 = and /3 = /q by doing a careful extrapo- 
lation to the thermodynamic limit, considering systems 
with up to 100 electrons. For (3 T-nP-Pur results agree 
with previous calculations of (5(iir).t30 
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FIG. 1. Energy gaps for creating neutral skyrmion pairs 
plotted vs. g. The number labelling each line segment is the 
total number of reversed spins S — 2K -\- 1. Results are given 
for thickness parameter (3 = and (3 — lo and have been 
extrapolated to the thermodynamic limit. 

When the Zeeman contribution is included the total 
energy gap for creating a neutral skyrmion pair for which 
the total number of reversed spins S = 2K -|- 1 is 



As{g)^A,f + 25{K) + Sg—. 



(7) 



These gaps are plotted vs. g in Fig. H for S' = 1, 3, 5 
and 7 and for /3 ^ and (3 = Iq. In contrast to Hartree- 
Fock calculationso which describe continuously varying 
spin textures, our results show the quantum nature of 
the skyrmion. As g decreases the energy gap undergoes 
level crossings in which the total number of reversed spins 
of the neutral skyrmion pair jumps by two. Note that al- 
though we have extended the S = 7 line down to g = 0.01 
we expect that as g decreases the density of level crossings 
will increase until the energy gap becomes, effectively, a 
smooth function of g. 

Increasing the thickness of the 2DEG is seen to lower 
the skyrmion energy gaps and favors excitations with 
fewer spin flips. This effect can be understood quali- 
tatively by noting that as the number of reversed spins 
K increases the charge of the skyrmion spreads out re- 
sulting in a lowering of its Coulomb energy. When the 
thickness is increased, and consequently the short-range 
part of the electron-electron repulsion is decreased, the 
excitation energies of the larger skyrmions will decrease 
less than those of the smaller skyrmions. 

We now turn to the effect of Landau level mixing 
on the skyrmion energy gaps. The importance of Lan- 
dau level mixing in a given system is characterized by 
the electron gas parameter r^ — l/{aB\/Tm), where n 
is the carrier density and ub = e?i /me^ is the effec- 
tive Bohr radius. It is straightforward to show that 
Ts — (ly /2y^'^ {e^ / elo) /huJc where hujc = heB/mc is the 



cyclotron energy so that, at fixed v, in the limit Ts ^ 
the cyclotron energy is much larger than the Coulomb 
energy and there is no Landau level mixing. However, in 
typical experiments Ts is of order 1 or higher and Landau 
level mixing cannot be ignored. 

To study Landau level mixing we Iw-ve used the fixed- 
phase diffusion Monte Carlo method,El which has proven 
to be a useful tool for studying 1|^ effect of Landau level 
mixing on quantum Hall states .ErlZl All of the wave func- 
tions considered here have been given in the form 

%p^ A V'spacc(ri, • • • , rw) «) (ii • • • Ik]k+i • • • Tjv) (8) 

where z/'spaco is completely antisymmetric under the ex- 
change of any pair of electrons i and j where either 
1 ^ «,i ^ -?^ or iC -|- 1 < i, j < TV. In a fixed-phase 
diffusion Monte Carlo simulation these wave functions 
are used as trial states by first writing the space part of 
the wave function in the form 

V'spacc(ri, • • • ,r^) = |^.pacc(ri, • • • , rAr)|e''^('-i'-''-~), (9) 

then 'fixing' the trial phase (^9, and finally constructing 
an effective bosonic Schrodinger equation for the positive 
definite wave function IV'spacol and solving that eqaation 
using standard diffusion Monte Carlo techniques.EJ The 
result of this procedure is the lowest energy state of the 
form (0) subject to the constraint that the phase is the 
same as that of the trial function. 



longer possible to transform a charge -\-e skyrmion into 
a charge — e skyrmion by particle- hole symmetry. There- 
fore, in order to calculate a physical quantity, we have 
computed the energy gaps for creating a charge ~\-e 
skyrmion with K spin flips and a single spin-flip quasi- 



electron. As before 
5{K) = EfP{K) ~ 
gaps are given by 



E^P and 
E^if{0), the corresponding energy 



letting As/ — 
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As{~g) = A,,f+SiK) + S~g—. 



(10) 



where now S = K + 1. 
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FIG. 2. Energy gaps for creating neutral excitations con- 
sisting of a charge +e single spin-flip quasielectron and a 
charge — e skyrmion with iC" = 0, 1 and 2 reversed spins for 
Vs — 0, 1 and 5 as a function of g for thickness parameter 
l3 — 0. Each line segment is labelled by the total number 
of reversed spins S = K + 1. Results are for a 30 electron 
system. 

We have implemented this procedure using ipgs, fpsf 
and ^sk{K) as trial wave functions to obtain, respec- 
tively, the fixed-phase energies Eg^ , Ef^ and E^if{K) 
for a system with 30 electrons. For finite r,, it is no 



FIG. 3. Spin and charge densities as a function of r/lo (r 
is the chord distance on the sphere from the center of the 
skyrmion) for a charge +e skyrmion with K — 1 for r^ = 
and Ts = 20. Results are for a 30 electron system. 

These energy gaps are plotted as a function of g in 
Fig. H for 5* = 1,2 and 3 and for r^ = 0,1 and 5. 
As Tg is increased the effect of Landau level mixing is 
seen to be qualitatively similar to that of finite thick- 
ness - the skyrmion excitation energies are lowered, with 
the smaller size skyrmions having their energies lowered 
the most. This effect can be understood qualitatively 
by examining the effect of Landau level mixing on the 
skyrmion spin and charge densities ps and pc defined by 



ps{r) 



{ijsk{K)\{p^{r)TP^{rmsk{K)) 



±p (11) 



where p^ is the number density operator for spin a and 
p is the uniform number density far from the skyrmion. 
Fig. |3| shows mixed estimatesEd of ps and pc for the K = 1 
charge -t-e skyrmion for rg — and rg = 20. For Ts = 20 
the possibility of mixing in higher Landau levels leads 
to a spreading out of the charge density of the skyrmion 
resulting in a lowering of its Coulomb energy at the cost 
of some kinetic energy. This effect is suppressed as K 
increases because for larger values of K the charge is al- 
ready well spread out and there is less energy to be gained 
by allowing the charge to spread further. 



In Fig. ^ the energy gaps are shown for 5* = 1, 2 and 3 
and r^ = 0, 1 and 5 for thickness parameter P = Iq. When 
finite thickness is included the Landau level mixing effect 
is seen to be much weaker than it is for zero thickness. 
For example, for g ~ 0.06, when /3 = the Ts = 1 energy 
gap is ^ 25% smaller than its rg = value, and when 
thickness is included the rg — energy gap drops ^ 45% 
for (3 = Iq, but the additional reduction of the energy gap 
due to Landau level mixing is only ~ 5% when r^ = 1. 
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FIG. 4. Energy gaps for creating neutral excitations con- 
sisting of a charge +e single spin-flip quasielectron and a 
charge — e skyrmion with K — 0,1 and 2 reversed spins for 
Vs — 0, 1 and 5 as a function of g for thickness parameter 
P = li). Each line segment is labelled by the total number 
of reversed spins S — K + 1. Results are for a 30 electron 
system. 



opposite is not the case. Ttiis is consistent with the ob- 
servations of Kralik et ahlij who performed variational 
Monte Carlo calculations of the effect of Landau level 
mixing and finite thickness on the v = 1 single spin-fiip 
energy gap A^^, finding that for rg = 1 the inclusion of 
finite thickness, of roughly the same size as that consid- 
ered here, significantly reduced the energy gap. Based on 
this observation Kralik et al. concluded that both finite 
thickness and Landau level mixing contribute equally to 
the reduction of the energy gap. However, we believe 
that if Kralik et al. had first included finite thickness 
and then Landau level mixing they would have observed 
the same strong suppression of the Landau level mixing 
effect reported here. 

To summarize, using skyrmion trial states recently in- 
troduced by MacDonald, Fertig, and Bre}^! we have per- 
formed variational and fixed-phase diffusion Monte Carlo 
calculations of the effect of finite thickness and Landau 
level mixing on skyrmion excitations at j/ = 1. We find 
that both effects lower the skyrmion excitation energies 
and stabilize skyrmions with fewer spin flips. However, 
we also find that the two effects do not work coherently 
together — when finite thickness is included the effect of 
Landau level mixing is strongly suppressed. 
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This weakening of the Landau level mixing effect due 
to finite thickness can be understood as follows. The 
thickness correction softens the short-range part of the 
Coulomb interaction, which in turn reduces the interac- 
tion energy the skyrmion stands to gain by delocalizing 
its charge. It follows that as /3 increases, the quasielec- 
tron charge delocalizes less for a given value of Ts, and 
so the reduction of the energy gap decreases. Because 
of this suppression of the Landau level mixing effect by 
finite thickness we beliesa; that the experimental observa- 
tion by Schmeller et al.Q of skyrmion energy gaps which 
are significantly smaller than theoretical estimates which 
ignore Landau level mixing is more likely to be due to 
disorder effects, which are poorly understood, than Lan- 
dau level mixing. Note that a similar reduction of the 
Landau level mixing effect due to finite thickness has 
been observed for spin-polarized Laughlin quasielectron 
and quasihole excitations in the fractional quantum Hall 
effect .0 

It is apparent from Figs. 2 and 4 that although finite 
thickness strongly suppresses Landau level mixing, the 
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